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Abstract 

It is critical to reduce the abnormal and complex jamming signals in the first order IMF component, because the decomposition 
result of the first order intrinsic mode function (IMF) in the empirical mode decomposition (EMD) method directly affects the 
subsequent decomposition effect, and also there are sampling restrictions for ultrasonic detectors. In this paper, we use a three- 
step optimization method to obtain the first order IMF. The first step is that the local extreme points of collected echo signals are 
fitted as the upper and lower envelopes by using the cubic spline interpolation method. The second step is to look for the curve 
of local mean of the echo signals on the basis of the upper and lower envelopes. The third step decomposes the curve of local 
mean by the EMD method. The simulation results show that our proposed method can effectively reduce jamming signals in 
the first order IMF and achieve the aim of improving the decomposition effect of the EMD method. 
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Introduction 

In modem industry, metal and nonmetal bonding materials are widely used in the National Defense Industry 
because of its excellent performance of composite. The quality of composite bonding materials has a decisive 
impact on production usage. One of many methods to identify the internal quality of products is the ultrasonic test 
technology. The principle of ultrasonic test is based on the propagation characteristics of ultrasonic in specimens. 
Firstly, the ultrasonic probe produces ultrasonic and uses a certain way to make it into specimen, then the 
ultrasonic travels along a straight line in the specimen. In the meantime it interacts with specimen material and its 
weaknesses. This makes the direction of propagation and characters of the ultrasonic changed. The ultrasonic 
whose direction of propagation is changed is received and displayed on the oscilloscope by the ultrasonic detecting 
device. This is called the echo signals. Then the echo signals are processed and analyzed by appropriate methods. 
Finally, the internal defects of bonding materials are evaluated according to the characteristics of the received 
ultrasonic. 

Echo signals detected by ultrasonic apparatus have non-stationary and nonlinear characteristics. In the face of such 
signals, they can be decomposed into a number of IMFs and a residual amount by using the EMD method. Every 
IMF contains the local features of signals in time domain and frequency domain, and the residual amount is 
expressed as variation trend of signals in time domain. In recent years, with the development of ultrasonic test 
technology, many new methods of signal processing are used for detecting structural damage. Since the Hilbert- 
Huang Transform (HHT) method based on the EMD has many advantages in dealing with nonlinear and non- 
stationary signals, this method has got the attention of researchers in the field of test engineering and got 
promotion application quickly. However, due to the EMD's own algorithm, the decomposition result of the first 
order IMF affects the subsequent decomposition effect. Therefore it is necessary to improve the decomposition 
result of the first order IMF. Taking into account that the EMD method is a data processing method based on the 
time domain, the accuracy of decomposition highly requires the data accuracy and smoothness. In this paper, an 
improved method based on the cubic spline interpolation is proposed, which represents a substantial improvement 
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over the original EMD. We use a three-step optimization method to obtain the first order IMF. The first step is that 
the local extreme points of collected echo signals are fitted as the upper and lower envelopes by using the cubic 
spline interpolation method. The second step is to look for the curve of local mean of the echo signals on the basis 
of the two envelopes. The third step decomposes the curve of local mean by the EMD method. The experimental 
analysis shows that the improved EMD method is very effective in the treatment of ultrasonic flaw detection. 


Principle of Empirical Mode Decomposition 


Empirical mode decomposition method, as a new method of digital signal analysis, was proposed by N.E Huang in 
1998, who was an academician of the National Aeronautics and Space in America. The EMD method decomposes 
signals according to the time scale characteristic of the data, which does not need to set any basis functions. 
Because of these characteristics, the EMD method can be applied to any type of signals in theory. In dealing with 
non-stationary and nonlinear signals, it has obvious advantages. The EMD is a "filtering method" to decompose a 
complex signal into a finite number of IMFs and a residual component. For a given signal, for example x{t) , the 
basic steps of the EMD method are as follows: 

(1) Identify all the local extreme points of signal x(t ) ; 


(2) Connect all these local extreme points with the cubic spline interpolation as the upper envelope u(t) and lower 
envelope v(t) ; 


(3) Obtain the curve of local mean of the upper and lower envelopes m{t ) = 


u(t) + v(t) 


(4) Calculate the deviation between x(t ) and m(t ) : h } (r) - x(t) - m(t ) ; 

(5) Treat h x (t) as a new original signal and repeat the above steps as many times as required until the envelopes are 
symmetric with respect to zero mean under certain criteria, then obtain the first order IMF component , that is 
c 1 (t) = h l (t) ; 


(6) The original signal x(t) minus c x (t) to continue repeating the above process until separate all IMFs, denoted by 
c x (t),c 2 (t),...,c n (t) , and list the following formula: 


x(t) = 'Lc l (t) + r n (t) 

i=i 

where r n (t) is the residue of signal x(t) and n is the total number of IMFs generated from the original signal x(t) . 
c i(t ) contains different frequency components of the original signal from high to low. 

Cubic Spline Interpolation 

Because the local extreme points of the original signals were fitted as the upper and lower envelopes by using the 
cubic spline interpolation in the EMD method, we use the cubic spline interpolation method to smooth the original 
signal in order to alleviate the error of fitting signal envelopes. The curve can better characterize the extrema of the 
original signals and achieve the purpose of improving the accuracy of EMD algorithm. 

The cubic spline interpolation is a smooth curve through a series of points. Contrasting to other methods, the 
advantage of the cubic spline interpolation method is to improve the smoothness of the piecewise linear 
interpolation function at nodes. For given different nodes x 0 ,x l ,x 2 ,...,x n , the corresponding values are 

Ten Ti > kh’--’ y n . pq e cubic spline interpolation function .S' (A) in the range of a =x 0 <.q <x 2 <...<. x n =b should 
meet the following conditions: 

(1) It has the second derivative in the range of [a ,b ]; 

(2) S(x k ) = y k (k=0,l,2,...,n) ; 

(3) S(x) is a thrice polynomial in each sub interval [x k ,x k+1 ] (k=0,l,2,...,n-l). 
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Simulation and Analysis 

In this paper, we use MATLAB programming software to construct a kind of cubic spline fitting method, which 
can smooth the echo signals. The reason to smooth the signals is that the collected signals are limited by the 
instrument which has the limitations of the sampling effect, then the amplitude of echo signals mutates more 
frequently (as shown in Fig. 1). If we use EMD method to process such echo signals directly, the first order IMF in 
the decomposition will have a lot of abnormal and complex jamming signals, which will affect the subsequent 
decomposition effect. Here we take an echo signal for example which comes from a thin board with a viscosity of 
30%. 
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FIG. 1 THE ECHO SIGNAL WITH A VISCOSITY OF 30% 
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FIG. 2 EACH ORDER IMF COMPONENT AND ITS SPECTRUM 
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Fig. 1 shows the echo waveform, which is collected by ultrasonic detector and does not do any processing. It can be 
seen from the figure that, the echo signal works as simple harmonic oscillation and it is made up of multiple 
echoes, which are gradually weakened with the increase of the propagation distance. The vertical coordinate of the 
echo signal stands for the amplitude of the wave, and it also represents the intensity of the echo signal that is 
produced by the ultrasonic which is emitted into specimen. The horizontal coordinate of the echo signal stands for 
the sampling points, and also represents the sound path. After decomposing the untreated echo signal by EMD 
method, the IMF components and the spectrum of each IMF are shown in Fig. 2. 

From Fig. 2 it can be seen that the frequency band of the first order IMF is wide, and there are many jamming 
signals in it. As the result of the first order IMF is the most influential to the following decomposition effect, and if 
we can reduce the error of the first order IMF, the decomposition accuracy of the whole data can be improved. So 
the local maximum points and local minimum points of the echo signal are fitted by cubic spline interpolation 
method, and upper envelope y (t) and lower envelope s 2 (t) are came into being. Then we obtain the curve of local 
mean of the upper and lower envelopes: 

40 = 

(as shown in Fig. 3, Fig. 4 shows the curve of local mean x(t ) ). Finally x(t) is decomposed by EMD method, and 
the decomposition process is shown in Fig. 5. 
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FIG. 3 SMOOTHING PROCESS OF THE ECHO SIGNAL 
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FIG. 4 THE CURVE OF LOCAL MEAN 
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FIG. 5 EACH ORDER IMF COMPONENT AND THEIR SPECTRUMS AFTER SMOOTHING 


66 


Signal Processing Research (SPR) Volume 4, 2015 


www.seipub.org/spr 


As seen from Fig. 5, the frequency band of the first order IMF is much narrower than it is shown in Fig. 2. That is to 
say, the number of jamming signals in the first order IMF is reduced after smoothing the signal by cubic spline 
interpolation method. The subsequent decomposition effect is better than it before. 


Evaluating Indicator 


In order to evaluate the accuracy of the IMF component between original EMD method and the new EMD method, 
in this paper, the correlation coefficient NC is chosen as the evaluating indicator. It is the degree of correlation 
between the IMF component and the corresponding signal. With this parameter the precision of IMF component is 
evaluated. Definition of the correlation coefficient NC is shown as (3-1): 


NC = 


N 


£ (x, - X)(y t - y) 

(=i 


ZCx,-xr»Z(y i -y) 


N N N 

N'Lx i y i -'£x i »'Ly i 

/=1 1=1 i=\ 

I N ~ N 7 | N " N 7 
>£*? -(£*,)• .NZyf-iZyf 
V 1=1 1=1 V i = 1 1=1 


(3-1) 


The performance comparison between traditional EMD method and the new EMD method is shown in Table 1. 
From the table we can see that, the NC of the first order IMF component obtained by the new EMD method is 
higher than that of the first order IMF component which is obtained by the traditional EMD method. That is to say, 
the accuracy of the first order IMF component obtained by the new EMD method is higher than that obtained by 
the traditional EMD method. In addition, the NC of the second and third order IMFs obtained by the new EMD 
method is lower than the corresponding IMFs obtained by the traditional EMD method. And the NC of the other 
IMFs are higher than those of the same order IMFs obtained by the traditional EMD method. It is further explained 
that the improved new EMD method is very effective in the ultrasonic testing technology. 


TABLE 1 COMPARISON OF PERFORMANCE EVALUATING INDICATOR OF DIFFERENT EMD METHODS 


correlation 
^^coefficient 
EMD methods^\^ 

NC(IMFl) 

NC (IMF2) 

NC(IMF3) 

NC(IMF4) 

NC (IMF5) 

NC (IMF6) 

NC (IMF7) 

NC (IMF8) 

NC (IMF9) 

Traditional EMD 

0. 6159 

0. 6811 

0. 2382 

0. 0169 

-0. 0031 

-0. 0005 

-0. 0126 

-0. 0098 

0. 0103 

Smoothing EMD 

0. 9203 

0. 2413 

0. 0465 

-0. 0213 

0. 0127 

0. 0334 

0. 0291 




Conclusions 

In ultrasonic test, the sampling effect of echo signals is limited by the instrument and other factors. The frequency 
band of the first order IMF obtained by EMD method is wide, and it comprises many jamming signals. In the 
traditional EMD method, the decomposition error of the first order IMF will directly affects the result of the 
following decomposition. In this paper, an improved method based on the cubic spline interpolation is proposed. 
We use a three-step optimization method to obtain the first order IMF. The first step is that the local extreme points 
of collected echo signals are fitted as the upper and lower envelopes by using the cubic spline interpolation 
method. The second step is to look for the curve of local mean of the echo signals on the basis of the two envelopes. 
The third step decomposes the curve of local mean by the EMD method. The experiment result shows that the 
optimization EMD method is significant for analyzing the breakdown signals containing complex and abnormal 
events. 
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